1988 


N8 9 - 1416 2 

* 7 

NA5A/ASEE SUMMER FACULTY RESEARCH FELLOWSHIP PROGRAM (j 


JOHN F. KENNEDY SPACE CENTER 
UNIVERSITY OF CENTRAL FLORIDA 


f 1 

ON THE HAZARD OF HYDROGEN EXPLOSIONS 
AT SPACE SHUTTLE LAUNCH PADS 



Prepared By: 

Academic Rank: 

University and Department: 

NASA/KSC : 

Division : 

Branch : 

NASA Counterpart: 


John M. Russell 

Associate Professor 

Florida Institute of Technology 
Department of Aerospace 
and Mechanical Engineering 

Mechanical Engineering 
Propellants and Gases 

W . 1 . Moore 


Date: 

Contract No.: 


August 19, 1988 


University of Central Florida 
NASA-NGT-60002 


207 



ON THE HAZARD OF HYDROGEN EXPLOSIONS AT 


SPACE SHUTTLE LAUNCH PADS 



by 

John M. Russell 
Florida Institute of Technology 


ABSTRACT 


This report was prepared in support of efforts to assess the hazard of acci- 
dental explosions of unburned hydrogen at space shuttle launch pads. It 
begins with a presentation of fundamental detonation theory and a review 
of relevent experiments. A scenario for a catastrophic explosion at a KSC 
launch pad and a list of necessary conditions fcontribuiing to it is proposed 
with a view to identifying those conditions which, if blocked, would prevent 
the catastrophe. The balance of the report is devoted to juxtaposition of 
reassuring and disquieting facts, presentation of a set of recommendations 
for further work, and listing of three main conclusions. One conclusion is 
that ignition of hydrogen-air mixtures by weak ignition sources in unconfined 
geometries may produce a detonation, provided the effective flame area in 
the initial fireball is rapidly increased by turbulent mixing. Another conclu- 
sion is that detonability limits can be different from and narrower than flam- 
mability limits only if one restricts the rate of work that can be done on 
a flammable gas by mechanical agencies acting on its boundaries. 


Key words: Deflagration-to-detonation transition, spherical explosions, 
detonation, hydrogen-air explosions, flame trench 
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SUMMARY 


The space shuttle main engines discharge hundreds of pounds of unburned 
hydrogen into the launch pad environment during shutdown from rated power 
level during an abort sequence and at the end of a flight readiness firing. 
Numerous personnel at KSC and elsewhere have conducted a variety of 
studies to determine whether this unburned hydrogen exposes the flight or 
launch hardware to the risk of damage from accidental explosions. Most of 
these studies were concluded or cancelled in late 1987. This report is t e 
result of a one-man ten-week effort undertaken in the summer of 1988 by 
a NASA/ASEE faculty research fellow. The main purpose of this effort was 
to make sense out of the results now in hand and to indicate areas in which 
future efforts might be especially valuable. 

In pursuing this goal. I found that few. if any, of the previous studies 
directed specifically to the problem of unburned hydrogen at space shuttle 
launch pads presented any kind of primer on the basic principles of mecha- 
nics. thermodynamics, and chemistry that govern the phenomena of detona- 
tion in gases. After some early discussions with my KSC contact (W.l. Moore 
of the Propellents and Cases Branch, DM-MED-4, KSC), we agreed that the 
preparation of a document which describes the basic ingredients of a detona- 
tion comments on the relevence of certain experiments on the origin of 
gaseous explosions to the conditions at KSC launch pads, presents definitions 
of all the relevent technical terms in related areas of combustion and deto- 
nation theory, and assembles a select bibiography would be worthwhile. 

This report begins with a definition of the problem and a description of the 
kinds of boundary conditions to which flow in the flame trenches at KSC 
shuttle launch pads are subject. After this, a connected account of the es- 
sential elements of the theory of gaseous detonations is presented. In this 
section, I review three older studies by G.l. Taylor on the theory of explo 
sions with spherical symmetry. As in the elementary theory of gasdynamics, 
much insight can be gained by considering the response of a gas to motion 
of a piston in contact with it, or to the analog of a piston in three dimen- 
sions (i.e. an expanding sphere). If a shock wave is sent through any 
flammable gas by movement of a piston or expanding sphere, then by moving 
the piston or expanding sphere fast enough, one may produce an arbitrarily 
large temperature rise across the shock and may thus cause the shocked 
gas to ignite. If, in addition, the motion of this piston or expanding sphere 
is sustained for a long enough time interval, then one can ensure complete 
combustion of the shocked gas. Such a zone of complete combustion following 
a lead shock amounts to a detonation wave. Overpressure produced by such 
detonation waves may easily reach tens of atmospheres, which illustrates their 
destructive potential. 

If a flammable mixture of gases is ignited at a point by a weak source, the 
thermal expansion of the resulting fireball may exert an action on the un- 
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burned gas around it similar to that of the expanding sphere described above. 
If the gas is a mixture of hydrogen and air, and if the flame front produced 
by weak ignition remains spherical, then according to results of a 1960 study 
by the Arthur D. Little Company the gas will not detonate. If, however, 
the effective area of the flame front is abruptly increased by turbulent mixing 
or by other mechanisms, then detonation may occur, as is indicated by several 
more recent experimental studies. 

After this discussion of fundamentals, the report proceeds to describe a scen- 
ario for a catastrophic explosion at a KSC shuttle launch pad and the neces- 
sary conditions contributing to it. The idea is that if a necessary contributing 
condition is blocked then the sequence of events would be interrupted and 
no catastrophe would occur. After a discussion of this catastrophe scenario, 

I juxtapose various reassuring facts with disquieting ones. The report then 
takes up the matter of recommentations and ends with a short list of conclu- 
sions. The recommendations are, first, to consider the list of necessary 
contributing conditions for a catastrophe discussed earlier and try to identify 
those that can be blocked most easily. The second recommendation is to 
exploit land, if necessary, develop) simple models of selected flow details 
to inform one's thought on the phenomena most likely to relate to these neces- 
sary conditions. 

There are three conclusions. In abbreviated form, they are, first, that deto- 
nations in unconfined hydrogen-air mixtures can be produced without blasting 
caps or other high energy detonators provided there is a mechanism for 
rapidly increasing the effective flame area in a local fireball. Secondly, deto- 
nability limits can be different from, and narrower than, flammability limits 
only if restrictions are imposed on the rate that work can be done on a flam- 
mable mixture of gases by adjacent boundaries. In the absence of such re- 
strictions, all flammabe mixtures are detonable. Thirdly, there are realistic 
prospects for ruling out accidental detonations of hydrogen at KSC shuttle 
launch pads. One promising approach is to establish a well-founded list 
of necessary contributing conditions for a catastrophe and to consider, 
through the use of simple models of selected flow details, the credibility or 
lack of credibility of each such condition. The results of such an effort 
would either certify the safety of existing equipment and operating proce- 
dures or identify changes that would lead to such a certification. 
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1. INTRODUCTION 


The main engines of the space shuttle orbiter discharge hundreds of pounds 
of unburned hydrogen into the launch pad environment during some modes 
of operation. Two such modes of operation are engine shutdown from rated 
power level in a flight readiness firing (FRF) and an abort. Millions of dol- 
lars have been spent by the U.S. Air Force and by NASA to assess the risk 
of a catastrophic accident resulting from explosions of this gas. Some of 
these efforts will be described in the following sections. As of July 1988, 
however, no one has demonstrated that clouds of explosive gases near the 
shuttle orbiter either definitely do or definitely do not constitute a credible 
hazard. 

In view of the importance of the unburned hydrogen problem, W.l. Moore 
of the Propellants and Cases Branch, Mechanical Engineering Division, Ken- 
nedy Space Center, NASA suggested to the present author that the prepara- 
tion of a general review of the unburned hydrogen problem would be a 
worthwhile project to work on during the author's NASA/ASEE Faculty Sum- 
mer Research Fellowship. This report is the result of this one-man ten-week 
effort . 

1.1 HISTORY 

According to internal NASA documents made available to the author, an acou- 
stic signature characteristic of a gaseous detonation was detected in the ex- 
haust duct of a space shuttle launch facility at the Vandenburg Launch Site 
in California sometime during (or before) November 1984. This event prompted 
a series of studies involving Rockwell International, Martin Marietta, Lock- 
heed Missies and Space Corporation, and technical personnel at several NASA 
and Air Force facilities (Vandenburg, Johnson Space Flight Center, Marshall 
Space Flight Center, the National Space Technology Laboratories, Kennedy 
Space Center, etc. ) which resulted in the design of a hydrogen disposal 
system for the exhaust duct at the Vandenburg facility [ cf. Breit & Elliott 
( 1987 ) 1 J The design featured a system to inject superheated water into the 
duct. The specifications call for a system in which the superheated water 
quickly turns to steam and is delivered in sufficient quantities to render 
the gases in the exhaust duct inert. Installation of this system was post- 
poned indefinitely after the Challenger disaster. 

The questions raised during the study of the Vandenburg facility raised sim- 
ilar questions regarding the safety of the shuttle launch pads at the Kennedy 
Space Center. Unlike the Vandenburg facility (from which no space shuttles 
have yet been launched) the KSC facilities have a distinguished performance 
record. In twenty five launches, six flight readiness firings and two aborts, 
no alarming pressure signatures have been reported The confidence one 
might derive from this observation is tempered, however, by the results of 
an elementary calculation due to F. Howard (1987) 2 at KSC. He noted that 
an explosion can damage the shuttle orbiter if it causes a shock wave to 


214 



strike the orbiter with an overpressure in excess of 1.32 psi. The expiosive 
yield of a spherical blast that produces this overpressure at the or b |te ' ' 
of course dependent upon the distance from the orbiter to the center 
the blast' Taking this distance to be 200 feet (which represents the distance 
from the base of the orbiter to the most remote part of the f ,am e trench a 
pads 39A and 39B at KSC) only six pounds of hydrogen are needed to fuel 
a blast that produces a damaging overpressure. 

During a shutdown from rated power level, by contrast, the amount .of un- 
burned hydrogen discharged by the space shuttle mam engine (after all 
hydrogen burn off ignitors have ceased to burn) is on the order of 
pounds. One may not, therefore, dismiss the problem out of hand ion the 
grounds that the amount of unburned hydrogen is too small to do harm. At 
this juncture, one should also remark that the figure 400 pounds refers only 
to the unburned hydrogen outside the space shuttle fuel tank. If detonation 
of that gas caused a breach of the shuttle fuel tank, the potential for harm 

could be far greater. 


Previous efforts to assess the hazard presented by unburned hydrogen have 
included : 


a. 


A statistical study of previous firings of hydrogen-oxygen rocket 
engines in the U.S. [Littlefield ( 1987) 3 ]. Among its results was the 
reassuring observation that only six confirmed detonations have taken 
place out of over 16,000 firings. 


b. Testing of a model of the exhaust duct at the Vandenburg Launch 
Facility into which exhaust gases from an H 2 -0 2 rocket were dis- 
charged in a controlled fashion. This work was done by the Lock- 
heed Missies and Space Corporation in Santa Cruz in support of the 
hydrogen disposal system for the Vandenburg Launch Site. 


c. Measurements of gas properties in the exhaust plume produced by 
firing a full-scale space shuttle main engine at the National Space 
Technology Laboratory in Mississippi. 


d Numerical simulation of the exhaust plume of a space shuttle main 
engine in various modes of operation. These calculations were per- 
formed by the Rocketdyne Corporation in September and October 

1986. 


e. Preparation to install instruments in the flame trench at pad 39B 

at KSC to measure gas properties during the flight readiness firing 
in the summer of 1988. 


f. Commissioning of a study by Bransford and Voth of the National Bu 
reau of Standards to assess the problem of unburned hydrogen at 
KSC and make recommendations [cf . Bransford & Voth (1987) J. 
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g. Consultation with personnel at Combustion and Explosives Research 
in Pittsburgh, Pennsylvania. This effort resulted in a recommenda- 
tion to install long-burning hydrogen burn-off ignitors in upcoming 
flight readiness firings. 

h. A study by E.E. Zukoski at the California Institute of Technology 
titled 'Flow into the SSME exhaust duct' (Aerospace Memorandum 
number 5 , February, 1987) and a reply by Shelby Kurzius (NASA 
MSFC) 

1.2 GUIDING QUESTIONS 

The author's efforts this summer have been guided by two questions, namely 

Q1 Is there a set of events common to all credible scenarios leading 
to a catastrophic explosion near the orbiter? Such a set of events 
might include: 

• formation of a suitably large cloud of explosive gas suitably close 
to the orbiter 

• initiation of a detonation in the cloud 

Q2 Is either of the above events precluded by present hardware and 
operating procedures at KSC? 

All other questions addressed in this report are stimulated by the two ques- 
tions just listed. 


1.3 OVERVIEW OF THE REMAINDER OF THIS REPORT 

In section 3, I will discuss the boundary conditions to which fluid in the ex- 
haust plume of the SSME is subject. Discussion of the physical phenomena 
that relate to gaseous detonations in the flame trench at KSC shuttle launch 
pads is not feasible without some background in the theory and observation 
of detonations. To this end, section 3 is devoted to an elementary discussion 
of the physical principles that govern detonation phenomena and to a dis- 
cussion of some of the most instructive examples from the theory of spherical 
explosions. Section 4 is devoted to a discussion of detonation experments. 

In section 5, I propose a scenario involving eleven events that lead to a cata- 
strophic explosion. Certain necessary conditions must contribute to these 
events and nine such conditions are listed in the same section. In section 
6, I juxtapose various reassuring and disquieting facts. Sections 7 and 
8 are devoted to recommendations and conclusions, respectively. 
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? boundary con djtions to which fluid in the exhaust plum e 

~ IS SUBJECT 

2.1 FIXED BOUNDARIES 

The fixed boundaries to which the exhaust plume from the SSME. is subject 
Jon^iJt of the mobile launch platform (with its exhaust duct), the blast de- 

these boundaries and their placement relative to each other. 

2.2 TIME-DEPENDENT BOUNDARIES 

The time-dependent boundaries to which the exhaust plume is subject consist 
If the hvdrogen burn off ignitors, the sound supression water spray (SSWS). 
the time dependent flow rates, chemistry, and thermodynamic state of t 
oases discharged by the SSME, and, in the case of a launch, the placement 
nf the SSME exhaust nozzles. Plots of the time-dependent cumulative dis 
Charge of unboned hydrogen from the SSME during -nous f, ring seguencea 
areoresented in Brcit 8 Elliott ( 1 987 ) 1 and Bransford 8 Voth ('9871 I to 
nama two sources I . The case of an on-pad abort appears to result m the 
greatest cumulative discharge of unburned hydrogen. Bransford and Vot 
„ iv _ thp fiaure 800 pounds of total discharge, though this figure aoes not 

are spent. The latter figure in the case of an on-pad abort 
of 400 pounds. 

None of the studies I have encountered presents reliable guantitadve informa- 
tion on the state of the gas beneath ( i.e downstream of) the SSWS. though 
some authors have stated conjectures about it. Dr. E.E. t 

writing as a consultant to the Aerospace Corporation prepared a docum 
titled 'Flow into the SSME Exhaust Duct 1 , Aerospace Me ^°^ an ^ U "\ "“'Tf er 
Febrarv 12 1987, in which he addresses this question. I have not y 

succeeded in obtaining a copy of this document, but I have read a review 
of it bv Dr Shelby Kurzius at MSFC. According to Zukoski and Kurzius, 
thJ^SWS discharges enough water so that, if thoroughly mixed with the ex- 
haust gases from the SSME, there will be enough to extinguish the exhaust 
flame by cooling. It remains an open question whether the exhaust flame 
is really extinguishable in practice and whether the downstream mixture of 
gases and water droplets is sufficiently d, luted with water vapor to render 

it nonflammable. 
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A - Volume of air-hydrogen sphere 
to produce 1.32 psig at SO*. 

B - Volume at 100'. 

C - Volume at 140‘. 

D - Volume at 200'. 



Figure 1. Space shuttle on the launch pad, viewed in a direction perpendi- 
cular to the flame trench [from Howard ( 1987) 7 ] showing relative placements 
of the exhaust nozzles, exhaust duct in the mobile launch platform, and the 
flame trench. Also shown are the results of some of Howard's calculations 
described in the text. 
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Figure 2. Space shuttle on the launch pad, viewed in a direction parallel 
to the flame trench. 
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3. BACKGROUND ON DETONATION THEORY 


As was stated in the Introduction, one can not discuss the specific issues 
relating to the hazard of detonations at KSC shuttle launch pads without 
having some working vocabulary on detonation phenomena. The present sec- 
tion furnishes two kinds of background that will be useful to this end, 
namely an analytical glossary of technical terms and a select catalog of some 
classical solved problems in gas dynamic and detonation theory that seem 
to be particularly informative. By 'analytical glossary', I mean a glossary 
of terms arranged in logical rather than alphabetical order starting with the 
concepts that are most fundamental. 

3.1 UNIDIRECTIONAL PROBLEMS 

3.1.1 INERT SHOCK WAVES. Consider a fluid flow in which the streamlines 
are parallel and the fluid speed does not vary in the cross-stream direction. 
Suppose, further, that the flow is steady. Let x be a spatial coordinate 
measured positive in the direction of the fluid velocity vector and let u be 
the component of the fluid velocity in this direction. Let p, p, and e de- 
note the pressure, mass density, and specific internal energy of the fluid. 

A shock wave is a discontinuity or step change in p, p, e, u, and other 
quantities that separates two regions in which these quantities are smoothly 
varying. The changes in p, p, e, u, etc. across a shock must be compatible 
with the basic laws of mechanics and thermodynamics including the law of 
conservation of mass, the equation for the rate of change of translational 
momentum, and the equation for the rate of change of energy. Let the sub- 
scrips '0' and '1' denote the conditions on the upstream and the downstream 
side of the shock, respectively. Then the three laws just mentioned may 
be expressed in the form 


e o + 


Po u o 

= Pl^l 

= m 

(1) 

Po + P 0 U 0 

= Pi + 

P 1 u 1 : 

(2) 

P o U D ? 


Pi , u i 2 

(3) 

Po + 2 

" e i + 

pi 2 


respectively. Here, m denotes the rate of transport of fluid mass per unit 
area (measured perpendicular to the streamlines). The parameters in these 
equations must be compatible with the thermodynamic equation of state of 
the substance. If the substance undergoes no chemical reactions in the 
shock itself, this equation may be represented by 

e = f(p.h 

on either side of the shock. If the substance is a thermally and calorically 
perfect gas, then we may write 
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in which y is a shorthand for the ratio c p /c v . Here, c p and c v are the 
specific heats of the substance at constant pressure and at constant volume, 
respectively. When equation (4) holds, equation (3) may be replaced by 


_1_ ?± + = _L_ El + iiil (5) 

y-1 Po 2 t -1 Pi 2 

The eariest discussion of shock waves appears to be due to George Gabriel 
Stokes [Stokes (1848) b ]. The foregoing equations are associated with the 
names of Rankine [Rankine (1870) 6 ] and Hugoniot [Hugoniot (1887) J and 
are often called the Rankine-Hugoniot shock conditions. 

In a typical problem in shock wave propagation, the value of y is a given 
constant. Also given are data specifying the thermodynamic state of the 
gas upstream of the shock, including the values of p 0 and p 0 . There are 
four unknowns, namely (u 0 , u x , p x . p x ) whose values are constrained by 
the three equations (1), (2), and (5). As one might then expect, exactly 
one parameter is left undetermined by the equations written so far. 

One convenient choice of parameter, in terms of which the family of solutions 
of ( 1 ) , (2), and (5) may be expressed is the Mach number based on the 
fluid speed u 0 and sound speed a 0 upstream of the shock, i.e. 



in which 



( 6 ) 

(7) 


The fractional changes in all the relevent flow quantities across the shock 
may then be expressed in terms of M, e.g. 

£l = 1 + -£<M’ - 1) (81 

Po ^r +1 

£i Hi ~ n+i ) m ? (9 ) 

P Q = Ul ‘ (y-DM 2 + 2 

as is shown in all elementary books on gas dynamics [cf . Liepmann £ Roshko 
( 1957) 8 ,§2 . 13] . One may regard M as a measure of the shock strength 
(stronger shocks having higher M). 

Now consider a plane which is oriented parallel to the shock and is situated 
downstream of it. Suppose, moreover, that this plane moves with the fluid. 
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Let x p (t) be the time-dependent position of this plane. Then the action 
exerted on the fluid in the region where x > x_ (t) (i.e. the region con- 
taining the shock) by the fluid in the region where x < x p (t) is indist- 
inguishable from the action of a rigid body such as a piston that drives 
the flow. If one imagines that the flow is in fact driven by a piston and 
that the motion of this piston is prescribed, then a suitable boundary condi- 
tion for the problem would assert that the velocity of the fluid in contact 
with the piston equals the (prescribed) velocity of the piston. Only when 
such boundary data are prescribed can one determine M (or any other mea- 
sure of the shock strength) uniquely. 

3.1.2 DETONATION WAVES. Let T denote the local instantaneous absolute 
temperature in a gas. From the equation of state of an ideal gas, i.e. 

p = ( c p“ c v) pT , 


it follows that 


Pi_ _ Pigi 
P o Pq^o 

Substituting (8) and (9) to eliminate p ; / p 0 arid Pj/pq , one obtains (after 
some rearrangement) 


T± _ Pi/Po 

T o Pi/Po 


2 ( T - 1 ) 

1 + n+u 2 


1M 2 +1 

M 2 


(M 2 -1 ) 


( 10 ) 


[cf. Liepmann S Rosko, equation 2.49]. If a gas is flammable and if a shock 
of sufficient strength passes through it, then either the abrupt rise in tem- 
perature, the abrupt rise in velocity, or both may initiate combustion of the 
gas. A definition of combustion will be furnished in section 3.1.3 below. 

If such a combustion takes place and is completed in a finite distance behind 
the shock, then the localized reaction zone and the shock that leads it is 
called a detonation wave. Detonation waves always propagate supersonically 
(i.e. M > 1) relative to the fluid upstream of them. A deflagration wave, 
by contrast, is a flame front that propagates subsonically relative to the 
unburned gas. 


Let the subscript '2' denote the conditions at the downstream extremity of 
the reaction zone of a detonation wave. The flow quantities upstream of 
the lead shock are related to those at the downstream end of the reaction 


zone by a set of equations similar to (1 

Po U 0 - 
P 0 + Po“o 2 = 


) - ( 3 ) , i.e. 


p ? u 2 = m 

(id 

P? + P2 U ? 2 

(12) 

P 2 U 2 

e ? + — + 4 

? p, 2 

(13) 
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Here, the specific internal energy e includes not only the usual thermal part 
(which for a calorically perfect gas equals c v r) but also a chemical part. 

If one assumes that the fluid on either side of the detonation is a thermally 
and calorically perfect gas and if one denotes by g > 0 the chemical energy 
per unit mas S Y that is converted to thermal energy during the chemical reac- 
tion one may take 

e n = (c v ) ( T 0 + g (14) 

_ / | T (15) 

e 2 ** ' cr v *::2 

in (13)- By manipulations similar to those which led to (5) above, equation 
(13) may be written 


Ho Po 
V 1 Po 


Y 2 P 2 t U 2 
^ 2 “ 1 p 2 2 


If values of ( Po . p„, Y 0 < Y a . q) are given and ( P? . p 2 , u 0 , u ) a ^e sought, 
then the system (11). (12). (16) consists of three equations for the four 
unknowns. As was the case with inert shocks, the family of solutions of 
this system will have one free parameter and again one may take this para- 
meter to be the Mach number M = u 0 /a 0 based on the fluid velocity an 

sound speed upstream of the lead shock. 


If one substitutes u 0 - m/p 0 ar, d u 7 
u c and Uj, one obtains 


m/p- into (12) and (16) to eliminate 


-4 - -U 

V P2 J 


P: ~ Po 

jf _L 

Po P 2 


^ ? P ? 

H."1 P 2 


Y o Po 
Yo _1 Po 


after some rearrangement. If one factors the expression ( p 0 . P* ^ j 11 
(18) and eliminates (m) 2 by means of (17). one obtains, after simplification. 


Y.-1 P? Ho" 1 Pi 


which involves only two of the four unknowns, namely p 7 and p 2 . A Hugo- 
niot curve is a graph of (19) with p 2 as the ordinate and 1/p 2 as the 
abscissa. This curve is a hyperbola [ cf . Strehlow (1968) , §5-4J. 

Now the speed of sound a 0 in the unburned gas is given by a 0 2 = YoPo/Po- 
It follows that 

( m) 2 = ( p 0 u 3 ) 2 = (p 0 a 0 M ) 2 = o 0 'Yo ~ = HoPoPqM 2 


2 2 3 



Thus, (17) can be written 


P2 

1 


- ^oPoPo M? 


( 20 ) 


- P 0 


±_ 

P: 


If p 0 , p c , Vj. and M are all specified, then a graph of (20) (with p 2 as 
the ordinate and 1/p 2 as the abscissa) is a straight line called the Rayleigh 
line . 

If a detonation wave is to exist, there must be at least one intersection be- 
tween the Rayleigh line and the Hugoniot curve. Depending on the strength 
of the lead shock (as characterized by the parameter M), there may be no 
intersection, one intersection, or two intersections of the Rayleigh line with 
the Hugoniot curve. 

I have already remarked that the strength of the lead shock (as characterized 
by the value of M) can be determined uniquely only if boundary conditions 
are specified ( e.g. by prescribing the motion of a piston that drives the 
burnt gas). The question of whether a detonation wave may or may not 
propagate through a given mixture of fuel and oxidizer can not therefore 
be settled without specifying some kind of boundary condition. 

If the boundary conditions are such that there is at least one intersection 
of the Rayleigh line with the Hugoniot curve, then the values of p 2 and 
p 2 at the intersection collectively define the thermodynamic state of the burnt 
gas. The overpressure associated with the detonation wave is defined to 
be p 2 - p 0 . 

The earliest attempt to calculate the speed of a detonation wave relative to 
the unburned gas by rational methods was apparently that of D.L. Chapman 
in England [ cf. Chapman ( 1899) 10 J. Chapman's results were independently 
rederived by E. Jouget in France [Jouget (1905,1906) 11 ]. Neither of these 
authors seemed to have fully appreciated the role of boundary conditions 
in determining the value of M in the system (19) and (20). Both authors 
instead proposed the ad hoc condition that the Rayleigh line intersect the 
Hugoniot curve at a single point of tangency. This condition is now called 
chapman- Jouget condition, or C-J condition. Determination of M by the C-J 
condition fixes the values of p 2 , p 2 , u 2 , etc. which are then called the 
chapman- Jouget values of these quantities. Now the C-J value of the speed 
of a detonation wave relative to the unburned gas is determined by the 
chemical and thermodynamic properties of the gases undergoing reaction. 

It thus furnishes a convenient reference scale in the presentation of experi- 
mental results. 

The first authors to substitute a proper treatment of boundary conditions 
for the C-J condition worked independently in three countries [ Y . B . Zeldovich 
( 1 940 ) 1 * in the USSR, J. von Neumann ( 1943 ) 1 3 in the U.S.. and W. Doring 
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( 1943 )- in Germany]. These authors also took the opportunity to incorporate 
finite reaction rate chemistry in modeling the structure of the reaction zone^ 

The contributions of these authors embody what is now called the Zeldovich 
von Neumann-Doring model of a unidirectional detonation wave. 

3 1 3 EXTERNAL WORK AND DETONABILI TY LIMITS. A chemical reaction 
involving oxidation may be slow or rapid. The prototype of a slow ox.dat.on 
reaction is corrosion (e.g. rusting of a sample of bare iron when exposed 
to air, moisture, and salt). Any common flame such as a candle flame may 
be taken as the prototype of a rapid oxidation reaction. The terms co us 
t ion and burning are synonyms that refer to any such rapid ox,dat '°": 
such rapid oxidation reactions are exothermic, i.e. they are accompan 
by a release of thermal energy. 

A homogeneous mixture of pure substances in thermodynamic equilibruim is 
called flammable if it will burn in response to some change in thermodynamic 
state (such as a rise in temperature). The flammability or nonflammability 
of a mixture of substances is, of course, dependent upon its chemical I compo- 
sition, i.e. the relative concentrations of the various substances in the m 
ture If a certain list of substances is specified and if one considers all 
the mixtures that may be formed from them, then some mixtures may be f am- 
mable and others may not. Any boundary that separates those mixtui es 
are flammable from those that are not is called a flammability limit . 

Flammability limits, as defined above, are chemical in nature and do not depend 
upon the boundary conditions to which a sample of material is subject. Some 
authors have proposed that the set of equilibrium mixtures that can be 
formed from a given list of substances can be divided unambiguously into two 
parts, namely 'detonable' and 'non-detonable' . If such a proposal wet e a 
good representation of nature, one could tabulate detonability limits from 
numerous detonation test results in the same manner as is done in the tabu- 
lation of flammability limits [ cf. Lewis S von Elbe (1961)-. table t 10 p 535]. 
The main difficulty with this proposal is that the changes in the f\u\d P™ 
perties across a detonation wave are not uniquiely determined by the chemica 
composition and thermodynamic state of the fluid approaching it. Indee , 
the strength of an ordinary (inert) shock wave, as may be characterized 
by the parameter M = u 0 /a 0 or the strength of the lead shock in a detonation 

wave is not determined entirely by the basic laws of mechanics, ^ her ^ n ° 
dynamics, or chemistry. The parameter M can. in fact, only be fixed by 
specifying boundary conditions, such as the motion of a hypothetical piston 
that abuts the fluid on the downstream side of the shock or detonation wave 
(as the case may be). 

Suppose that a flammable mixture of gases will ignite if its temperature is 
raised above some threshold value T ai (i.e. the autoignition temperature |. 

Bv sending a piston-driven shock wave through it and by prescribing a pis- 
ton speed sufficient to raise the ratio T /T 0 [cf . equation (10) above] above 
the level r ai /T 0 , one may ensure that the gas will ignite and that the 
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resulting reaction zone will follow the shock. Such a procedure amounts 
to producing a detonation wave and the above argument suggests that all 
flammable mixtures are detonable. 

The argument given above hinges on the idea that a piston can do work on 
a body of fluid in contact with it. If one allows the piston to do an arbitra- 
rily large amount of such work in a given time interval, then one can always 
send a detonation wave through any flammable mixture of gases. This argu- 
ment implies that 'detonability limits' are as broad as 'flammability limits'. 

If. however, one incorporates restrictions on boundary conditions into one's 
definition of detonability limits, then one may well arrive at detonability limits 
that are narrower than flammability limits. One such restriction might assert 
that the piston is at rest relative to the unburned gas far ahead of it. 

3.2 MULTIDIRECTIONAL PROBLEMS 

3.2.1 INERT SHOCK WAVES 

3. 2. 1.1 Taylor's expanding sphere problem . As was pointed out in section 
3.1.1, the one-parameter family of solutions of the equations for steady non- 
reacting unidirectional flow with a shock wave can be generated by the solu- 
tion of a 'piston problem'. Specifically, if gas in a tube is initially at rest 
and if the gas is set in motion at t - 0 by a piston whose velocity rises 
abruptly in the manner of a step function, then a shock wave propagates 
ahead of the piston into the region of undisturbed gas. The strength of 
this shock can be represented by the Mach number M = u 0 /a 0 (see section 

3.1.1 for definitions of the symbols). This shock strength, in turn, can 
have any positive value. Its value is determined by the speed of the piston. 
The values of all of the flow quantities behind the shock are then determined 
by equations (8), (9), and (10). 

In the late nineteen thirties. C.l. Taylor in England addressed the question 
of whether one can formulate a problem in spherical geometry that is analo- 
gous to the piston problem described above. Taylor's efforts resulted in 
a manuscript titled 'The air wave surrounding an expanding sphere' which 
was submitted for publication to the Royal Society of London in 1939. It 
was not published until 1946 [Taylor ( 1 946 ) 1 6 ] . owing, apparently, to secu- 
rity classification during World War Two. 

Taylor's results show that if the radius R of the expanding sphere increases 
in time at a constant rate u 2 . then a shock wave propagates into the sur- 
rounding fluid. The radius r t of this spherical shock also increases with 
time at a constant rate. T aylor's results include tabulations of the distri- 
butions of u/a and p/p 0 versus r/(at ) and p = U 2 /a , in which a and p 0 
are the speed of sound and the fluid pressure, respectivley, in the remote 
undisturbed air, u is the local instantaneous fluid speed in the region be- 
tween the expanding sphere and the shock, r is the distance from the cen- 
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ter of the expanding sphere to any place where values of the flow quantities 
are sought, t is the time (relative to a hypothetical reference time when 
the expanding sphere had zero radius), and u 2 (as stated earlier ) is the 
velocity of expansion dR/dt of the expanding sphere. As before, let the sub 
script M' denote conditions immediately downstream of the shock. The ratios 
Pl / Po and r 1 /r 0 may be expressed in terms of the ratio p x /p 0 by means o 

the formulas 

Pi _ pr-n + (t+i ) (Pi/Po ) d £1 . Pi 

p 0 ‘ 0+1 ) + O-lMpi/Po) T o P o Pi 

the first of which is deduced from the system (1), (2), and (5) [ibid.. 
equation (24)] and the second of which follows from the equation of state 
of an ideal gas. Table 1 and Figures 3 and 4 illustrate some of Taylors 
results. In each of the smaller tables in Table 1. the value of p/po corres- 
ponding to the largest value of r/(a t) is equal to Pl /p 0; Thus, I aylor s 
tables permit determination of Pl , Pl . and ^ in terms of the giving para- 
meter p and the thermodynamic state of the undisturbed gas. The distr b 
tions of p/p, and t/t x between the shock and the expanding sphere then 
follow from the tabulated values of p/p 0 and the adiabatic relations 

(-Y-l)h 


p/Pi = (p/pj 


1/y 


and 


T/T j = (p/pj 


These results of Taylor would seem to offer promise in modeling the imtia 
tion of detonations in spherical geometry. Consider, for example, a localized 
region of burning gas in which the flame surface (or surfaces) are deflagra- 
tions. The instantaneous rate at which thermal energy is released by the 
burning is influenced by many variables, one of which is the instantaneous 
area of the flame front. For the purpose of comparison, suppose that the 
flame front at some instant is a closed surface that encloses a given volume. 

A flame front in the shape of a smooth sphere would have a smaller area 
(and would, therefore, yield a smaller rate of energy release by burning) 
than would a flame front with any other shape. At the opposite extreme, 
one can imagine a flame front with a highly convoluted shape (perhaps as 
a result of turbulent mixing) whose rate of energy release is, say, two 
orders of magnitude higher than is that of the spherical flame. Such a region 
of burning gas would, of course, expand and the flow about Taylor's 
expanding sphere might be a suitable model to describe the action of a local- 
ized region of burning gas on the unburned gas surrounding it. If one 
supposes that such a model is appropriate, then a scenario for the initiation 
of detonation might run as follows. The flame front in a localized region 
of burning gas becomes highly convoluted as a result of turbulent mixing. 
This convoluted flame front results in a release of thermal energy that is 
much larger than what would be released if the flame front were smooth and 
spherical. This rapid release of thermal energy causes the ball of burning 
gas to expand. This expanding ball exerts an action on the unburned gas 
surrounding it in the manner of Taylor's expanding sphere including, in 
particular, the production of a spherical shock wave. The temperature and 
velocity of the unburned gas both suffer abrupt rises as the spherical shock 
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Table 1. Results from Taylor (1946) 16 on the air-wave surrounding an expan- 
ding sphere. Here a u : /c 2 , in which c 2 is the speed of sound in the 
air at the surface of the sphere. See text for definitions of the remaining 
symbols. Note that p = u 2 /a. in which a is the speed of sound in the remote 
undisturbed air. Thus, a and p are two alternative parameters for represen- 
ting the rate of driving of the air by the expanding sphere. 
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Figure 3. Plots of the pressure distributions in Table 1 [Taylor ( 1 946 ) 1 6 ] . 
Each curve corresponds to a particular rate of expansion of the sphere. The 
leftmost point on each curve corresponds to the surface of the sphere; the 
discontinuity at the right of each curve corresponds to the shock wave. 
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Figure 4. Plots of the radial velocity distributions in Table 1 [Taylor 
(1946) 16 ]. As in Figure 3, the step discontinuities represent the shock wave 
for each rate of driving. The dashed line represents the trajectory of the 
surface of the expanding sphere. 
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oasses by If the latter is of sufficient strength to satisfy the conditions 
for ignition of the shocked gas, then the shocked gas will begin to react 
If moreover, those conditions are sustained for a sufficient duration < D f t 
then the reaction will proceed to completion. Such a shock followed by a 

reaction zone is. of course, a detonation wave. . T |"; ^asserted 
constitutes a mechanism for initiation of a spherical detonation, as asserted. 

3. 2. 1.2 Taylor's fixed-total- energy blast wave problem. In a review art. - 
He published in Annual Reviews of Physical Chemistry , J.H.S. Lee or McGill 
University [Lee (1977) 17 ] discussed efforts to model the intiation of spheri 
cal detonations. One such effort (cf. p 91 of that paper) «P^ e , d Lher 
of C.l. Taylor's contributions to the theory of explosions [T a y ,o r ( 1950) J. 
This Daoerwas written in 1941 in support of the early work on the atomic 
bomb ^and remained classified for nine years. Taylor formulated the problem 
by supposing that a spherical blast wave is generated by the sudden release 
of a fixed total energy E , which then represents the sum of the thermal and 
kinetk energies of thl blist. In contrast to the action of conventional high 
explosives this release of energy was not accompanied by a release of gas 
through the vaporization of condensed matter. As in the expanding sphere 
problem, Taylor sought a solution of the equations of motion of compressible 
Unviscid nonheatconducting) flow in which the distributions of the *'° w 
Quantities expressed in terms of appropriate nondimensiona variables, are 
self- si m Mar at all times. Taylor showed that such self-similar solutions exist 
only in the case when the strength of the shock wave is asymptotically large, 
thus permitting one to replace the rightmost member of (9), for example, 
by its asymptotic limit as M -*■ « . 

T aylor's results include formulas relating the radial velocity JR/dt of _ the 
shock to the time t for a given blast energy E. In the model discussed 
by Lee a critical value of E was proposed by equating the Chapman- Jouget 
velodty of a detonation wave [ cf. section 3.1.2 above] with dR/dt and equa- 
ting the time t to the time t for completion of the chemical reaction (the 
so-called induction time). Lee remarks that the values of E ^ compu ^ 
are about three orders of magnitude less than are the energies of b lasting 
devices needed to produce spherical detonations in experiments. Lee off* ers 
some ad hoc explanations for the discrepancy. Curiously, Lee does " ot 
impute significance to the discrepancy between the necessarily moderate shoe 
strength in the shock waves produced by blasting caps and the infinite i shock 
strength assumed in Taylor's fixed- E blast wave model. As Taylor points 
out in his paper, his fixed- E blast wave solution is comparable to blasts p 
duced by conventional condensed matter explosives only if the mass of the 
air enclosed by the shock is much greater than the mass of the explosive. 

I do not believe that the conditions of the blasting cap experiments cited 
by Lee fulfill this condition within the time interval of interest any better 
than they fulfill the condition of asymptotically large strength. 

Lee (1977) does not mention any of Taylor's papers other than the one con- 
taining Taylor's theoretical model of the atomic bomb blast. He may be aware 
of Taylor's work on the air wave surrounding an expanding sphei e and have 
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good reasons for rejecting it as a model for initiating spherical detonations. 

I have not seen any such reasons, however, and lacking them, I am inclined 
to regard the mechanism for initiating spherical detonations described in the 
preceeding section as more believable than the one described by Lee (1977). 

3.2.2 PROPAGATION OF A SPHERICAL DETONATION WAVE. In his many 
studies of spherical explosions during World War Two, Taylor also addressed 
the problem of how to model the propagation of a spherical detonation wave. 
Taylor's theory of the spherical detonation was formulated in the same year 
as and prior to his work on the fixed total energy blast wave. At that time 
(January 1941), the standard model for the propagation of detonation waves 
was the Chapman-Jouget model discussed in section 1.2 above. As I stated 
there, the important role of boundary conditions in determining the shock 
strength is not incorporated in the C-J model. Taylor's spherical detonation 
wave is predicated on the assumption that the C-J condition is satisfied and 
the applicability of Taylor's spherical detonation model is limited accordingly. 
In their book Detonation, Fickett and Davis ( 1 979 ) 2 0 remark that the problem 
of spherical detonations has still not been properly treated. It is possible, 
therefore, that Taylor's World War Two contribution, limited as it is, had 
not been superceded as of 1979. At this point, however, one should call 
attention to the important book Similarity and Dimensional Methods in Mechan- 
ics by L.l. Sedov (Sedov (1959) 21 J. In chapter four of that work, Sedov 
formulates a general analytical framework for generating self-similar solutions 
of the equations of gas dynamics in one, two and three dimensions. Sedov 
not only recovers all of the results of Taylor that I have cited so far, but 
is also able to replace some of Taylor's numerical solutions with closed-form 
analytical ones. Sedov's contribution goes beyond Taylor's in that Sedov is 
able to delineate the complete set of circumstances under which self- similar 
solutions of the equations of reactive gas dynamics with one space coordinate 
(as occur, for example, in the problem of spherical detonations) are pos- 
sible. Fickett and Davis do not cite Sedov's book, so it may be that their 
characterization of post war work on spherical detonations is overly harsh. 

Taylor's model of the spherical detonation wave is consistent with the view 
that a spherical detonation, once initiated, can propagate through the whole 
region occupied by explosive, i.e. there is no a-priori reason to suppose 
that a spherical detonation will extinguish itself after it reaches a certain 
radius. 

As in problems discussed earlier, Taylor sought and found a solution of the 
equations of motion of a gas in which the distributions of the flow quantities 
expressed in terms of appropriate nondimensional variables, were self-similar 
for all times. Since the details of these distributions are less germane to 
the purposes of this report than is the fact that the whole cloud of gas deto- 
nates, I will set aside further discussion of Taylor's theory of spherical deto- 
nations. 
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3.3 REMARKS ON CONFINEMENT 


The possibility that a detonation wave will ultimately result from ignition 
of a given sample of flammable gas is strongly dependent upon boundary 
conditions. Thus, if the gas is in a tube closed at one end and if the gas 
is ignited at the closed end, then the likelihood that the resulting flame front 
will evolve into a detonation is much greater that if the gas were in a spheri 
cal balloon high above the ground and the gas were ignited at the center. 
These two geometries typify 'confined 1 flows and 'unconfined flows, respec 

tively . 


The notion of confinement is hard to quantify directly, at least if one tries 
to tie it to geometries of particular solid boundaries in contact with an explo- 
sive gas. What really seems to matter is the dimension of the space in which 
the gas is allowed to move. Thus, if a gas is constrained to move along 
parallel streamlines, it is more confined that if it is allowed to move radial y 
along rays perpendicular to an axis ( i.e. in two dimensional cylindncally 
symmetric motion). This motion, in turn, is more confined than is motion 
along rays emanating from a point ( i.e. in three dimensional spherically 
symmetric motion). 


The idea that one can better gauge confinement by counting space dimensions 
in which gas may move than by looking at the detailed geometries of confin- 
ing walls becomes clear when one contrasts the case when gas is ignited in 
a tube closed at one end with the case when gas is ignited in a tube closed 
at both ends. If confinement is to be implicated as a factor that always 
increases the likelyhood of detonation and if wall geometry were the essence 
of confinement, then gas in a tube closed at both ends should be more prone 
to detonation than is gas in a tube closed at only one end. I do not believe 
that such a prediction would be borne out be experiment. Venting of the 
tube at one end allows the burning gas at the other end to expand and act 
like a piston that sends a shock wave ahead of it. The formation of this 
shock wave is a basic step in the evolution of a detonation wave. Sealing 
a tube at both ends could inhibit turbulent mixing of a initial subsonic flame 
front, thereby preventing it from accelerating to supersonic velocity and 
thus inhibiting one mechanism for shock formation. 


tK BACKGROUND ON DETONATION OBSERVATIONS 

The preceeding section was devoted to descriptions of basic physical pheno- 
mena, definitions of terms, and attempts to isolate the basic cause and effect 
relationships in such phenomena as the initiation of unidirectional and spheri- 
cal detonation. Thus, while the last section was devoted primarily to theory, 
the present one is devoted to observations. The number of new papers ap- 
pearing each year in which results of experiments on detonation are reported 
is quite large. I have found the review articles by J.H.S. Lee 1*9/7) 
and Marshall Berman ( 1985) 22 to be especially useful as introductions to 
this vast literature. The book by Strehlow (1968) 9 is also an excellent 
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introduction to the general subject of combustion theory and observation. 

The experimental studies discussed in this section are restricted to those 
I have read that seem to me to be most informative in regard to the hazard 
of hydrogen explosions at KSC shuttle launch pads. 

4.1 PROPAGATION OF DETONATION WAVES DOWN TUBES AND CHANNELS 

4.1.1 LOCAL EXPLOSIONS. Consider a tube filled with an explosive mixture 
of gases and closed at one end. It may happen that if the gas in the tube 

is ignited at the closed end. a flame propagates toward the open end, aceler- 
ates and eventually evolves into a detonation wave that propagates all the 
way to the open end of the tube. A beautiful set of photographs, which 
record many of the physical processes that take place in such an experiment 
was published by Urtiew and Oppenheim ( 1966) 23 . One curious feature of 
detonation waves captured by these photographs is the repeated occurence 
of localized spherical explosions (which Urtiew & Oppenheim called 'the explo- 
sion within the explosion'). Thus, even though the lead shock in a detonation 
wave may be very nearly planar, the flow behind it is, by no means unidi- 
rectional. The intermittent generation of three dimensional local explosions 
seems to accompany all observations of detonation initiation and propagation. 

4.1.2 DETONATION CELL WIDTH AND THE EMPIRICAL FORECASTING OF 
DEFLAGRATION-TO-DETONATION TRANSITION. If the inside surface of 
a tube is covered with a smoked foil and if a detonation wave is then sent 
down the tube in the manner described above, then the detonation wave will 
scour the foil selectively leaving a cell-like pattern on it [of. Figure 5, taken 
from Strehlow (1968) 9 J. The task of determining the precise mechanism 

to account for this wall scouring is a challenge for theoreticians that need 
not concern us here. Whatever their cause, however, the detonation cells 
furnish an experimentally observable length scale that typifies the thickness 
of a detonation wave. 

The detonation cell width depends upon the mole fractions of the fuel, the 
oxidizer, and of any diluting substance that may be present in the mixture 
of gases upstream of the detonation wave. Given a particular chemical 
reaction involving combination of a fuel with an oxidizer, one may define 
a special ratio of the molar concentration of fuel to the molar concentration 
of the oxidizer which corresponds to the ratio of the same constituents in 
the reaction product. Such a ratio is called stochiometric . Thus, in the 
reaction 2H 2 + o 2 -* 2H 2 , a stochiometric mixture of h 2 to o 2 would contain 
two moles of h 2 for each mole of o 3 . A parameter that represents the close- 
ness of a given mixture to a stochiometric one is the equivalence ratio 4> 
defined as follows: 

mole fraction of fuel to oxidizer in a "N 

given mixture of substances ]_ 

mole fraction of fuel to oxidizer in a stochimetric ^ 
mixture formed from the same substances , 


'equivalence^ 
^ ratio J 


234 


qTUGINM' 

ok ?^° K 


‘ ^ lTY 



Figure 5. Photograph of the pattern produced on a smoked foil on the wall 
of a tube down which a nominally planar detonation wave propagates [from 
Strehlow (1968) 9 ]. The wave propagated from left to right. The symbol 
X denotes the detonation cell width. 
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Thus, a mixture is fuel-rich or fuel-lean according to whether $ is greater 

than one or less than one. respectively. The value $ = 1 corresponds, of 
course, to a stochiometric mixture. 

By plotting detonation cell width versus equivalence ratio, one has a basis 
for comparing results involving different chemical reactions. An example 
of such a plot is given in Figure 6 from Berman (1985) 22 , which illustrates 
the dependence of detonation cell width (measured parallel to the plane of 
the detonation wave) upon equivalence ratio $ for hydrogen-air detonations 
with varying degrees of steam dilution. There does not appear to be any 
simple accurate fundamental method for the direct analytical calculation of 
the detonation cell width. By allowing for finite reaction rates in a one di- 
mensional model of the reaction zone behind the lead shock, however, one 
can define a length scale (the induction length) that scales with the thick- 
ness of this zone. This induction length may be correlated with the detona- 
tion cell width and the curves in Figure 6 illustrate such a correlation. 

Many empirical correlations may be expressed in terms of the detonation cell 
size. Thus, for example, if the plotted detonation cell size for a given mix- 
ture of gases is larger than the bore of a tube in which one wishes to deto- 
nate that gas, then there is reason to believe that in that tube the gas (in 
the absence of piston driving) will seem less detonable. Several authors 
have proposed such ideas and have reported evidence in support of them. 

One remark made by many such authors is that detonability limits are scale- 
dependent and that larger-scale boundaries are more conducive to detona- 
tion than are boundaries of smaller scale [ cf . Berman ( 1 985 ) 2 2 ] . 

A second kind of empirical correlation involving the detonation cell width 
is a correlation between the size of the smallest obstacle capable of influen- 
cing deflagration-to-detonation transition (DDT) for a given mixture of gases 
and the standard width of the detonation cells for that mixture. Sherman, 
Tiezen, and Benedick (1986) 2t * , for example, have found experimental evidence 
that a deflagration front propagating down a channel across a regular array 
of fence-like obstacles may undergo transition to detonation and that this 
transition to detonation will be influenced by the fences if the spacing 
between them is of the order of two to three times the detonation cell width. 
The validity of such a correlation would indicate that small detonation cell 
sizes correspond to more explosive gases, i.e. gases in which DDT may 
be triggered by smaller obstacles. For stochiometric mixtures of hydrogen 
and air (with no steam dilution) the detonation cell width A is between four 
and six millimeters [ cf. Figure 6]. If, however, the w 2 -o 2 ratio in the 
mixture is kept the same but steam is added to bring the total steam fraction 
to twenty percent, then the detonation cell width rises to a value between 
ten and twelve centimeters. Raising the steam fraction to thirty percent 
raises the detonation cell size to about thirty centimeters. Hydrogen-air- 
steam mixtures with steam fractions above fourty-five percent are so inert 
that no attempts to produce detonation waves in them have yet succeeded. 

The observations on detonation described above suggest that the addition 
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Figure 6. Detonation cell width X versus equivalence ratio <J> for hydrogen- 
air-steam mixtures and correlation with calculated 'induction lentgh' (i.e. a 
typical streamwise thickness of a planar detonation wave, as may be calcu- 
lated in the Zeldovich-von Neumann-Doring model) [from Berman (1985) 22 ]. 
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of steam to the exhaust plume of the space shuttle main engines exerts a 
powerful effect in reducing the detonability of the gases in it. Such consi- 
derations undoubtably played a role in the design of the steam inerting 
system for the exhaust duct at the Vandenburg launch site. In this context, 
the effect of the sound supression water spray (SSWS) on the main engine 
exhaust plume is obviously favorable. 

4.2 EXPERIMENTS ON UNCONFINED DETONATIONS. An older study by 
the Arthur D. Little Company [Anonymous (I960) 25 , hereinafter referred 
to as ADL] contains an impressive compendium of engineering work done in 
support of the Atlas Centaur program. The report contains results of exper- 
iments on detonation of hydrogen-oxygen-air mixtures in spherical balloons 
situated at a remote distance [ i.e. many balloon diameters] away from the 
ground or any other walls. In these experiments, detonations were produced 
by placing some flame source at the center of the balloon and causing it to 
ignite. The flame sources ranged from 'weak* sources, represented by 
sparks, to 'intermediate' sources, represented by squib flames, and 'strong' 
sources represented by blasting caps and other condensed matter detonators. 

Some of the results reported in ADL seem reassuring in the context of this 
report. Thus, the authors of ADL report (p 18) that 'detonation will 
not occur when hydrogen is vented to the atmosphere as long as the hydro- 
gen-air mixtures are unconfined and are initiated by ignitors of the non- 
explosive type'. This conclusion was based upon observations of spherical 
flames produced at the centers of balloons by weak sources. In these exper- 
ments, no direct means were provided to produce convolutions of the flame 
fronts and it may well be that the flames remained nearly spherical. As was 
remarked in section 3. 2. 1.1 above, a spherical flame front has the smallest 
surface area and concomitant rate of energy release of all possible flames 
that enclose a given volume. Thus, the relatively benign characteristics 
of the spherical flame produced in the ADL experiments on hydrogen-air 
mixtures may not be representative of what could happen if the flame front 
were allowed to be wound up in turbulent eddies (as might be produced, 
for example, by passage of the flame front across an array of turbulence- 
producing obstacles). The authors of the ADL report were aware of the 
importance of turbulence in the flame front. Indeed, they cited such turbu- 
lence as their reason for abandoning an attempt to determine the minimum 
spark energy necessary to trigger a detonation [cf. ADL, pp 57-58]. One 
may surmise from the discussion in ADL that the direct effects of the spark 
were overwhelmed by the effects of turbulence produced by the instrument 
supports. 

4.3 FURTHER OBSERVATIONS ON DEFLAGRATION-TO-DETONATION TRAN- 
SITION. The article by Berman ( 1 985 ) 2 2 reviews several experiments on 
the propagation of flames through bags filled with flammable gases. In one 
example [ cf . Moen, Bjerkvedt, Jenssen, and Thibault (1985) 26 ], a flame 
front propagates down the bag, undergoes distortion, and forms a 'tongue 

of flame' that advances ahead of the main flame front ( cf Figure 7). When 
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this tongue of flame reaches the far end of the bag, it is wrapped up in 
a corner eddy. A detonation then originates within this corner eddy and 
the resulting shock wave propagates through the rest of the facility. 

In this experiment, the gas bag abutted a solid floor, but was otherwise 
unconfined. Although the mixture of gases was acetylene and air, Berman 
(1985) 2 2 points out that the mixture has the same equivalence ratio and is 
thus comparable to a hydrogen-air mixture with 19% hydrogen. For compari- 
son, a stochiometric mixture of hydrogen and air has 30% hydrogen. The 
detonation observed by Moen, et al. was not planned as part of their exper- 
iment. The reassuring conclusion from the ADL report cited earlier loses 
some of its force when compared with these results of Moen et al . 

Results analogous to those of Moen et al . have been obtained by German 
investigators [see Berman (1985) for the original sources] in a large gas 
bag supported by a metal frame. It appears that passage of a deflagration 
over one of these supports resulted in the shedding of an eddy which then 
became the center of a spherical detonation. 

In assessing the relevence of these experiments to the conditions at KSC 
shuttle launch pads, it is important to realize that the gases in the gas-bag 
experiments were well-mixed, the bags prevented the gases from dispersing, 
and the gas was initially at rest. All subsequent motion was due to non-uni- 
form energy release by burning. At the KSC shuttle launch pads, the gases 
are not well-mixed, they are free to disperse, and there is considerable mo- 
tion in the flame trench, owing to the momentum of the gas in the deflected 
SSME exhaust plume. The flow inside a jet engine combustor might, in fact, 
represent the conditions in the flame trench more accurately that does the 
flow in the gas-bag experiments, owing to the nonuniformity of the fuel- 
oxidizer mixture, the presence of an overall through flow, and the presence 
of turbulent mixing. 

SCENARIO FOR A CATASTROPHIC EXPLOSION 
5.1 A SEQUENCE OF ESSENTIAL EVENTS 

In section 1.2 above, I listed two questions that guided my work this summer, 
the first of which was 'Is there a set of events common to all credible 
scenarios leading to a catastrophic explosion near the orbiter? 1 For the pur- 
pose of discussion, I will postuate that the answer is 'yes' and propose the 
following sequence of essential events. 

El One or more of the space shuttle main engines discharges unburned 
hydrogen . 

E2 The usual flame in the exhaust plume is extinguished (at least 
locally) . 

E3 Some of the hydrogen in the exhaust plume (perhaps the gas below 
the water curtain) collects in an unreacted cloud rather than burning 
in the usual blow-torch fashion. 
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E4 At some later time, the gas in the unreacted cloud is ignited. 

E5 The flame front is subjected to rapid distortion by turbulent mixing, 
thus causing a disproportionately large rate of release of thermal 
energy compared to that released by a laminar flame. 

E6 Thermal expansion of the rapidly burning gas forces the unburned 
gas in its neighborhood out of its way, thereby sending a compres 
sion wave into it. 

E7 The compression wave generated in E6 steepens to become a shock 
wave in the usual way. 

E8 Passage of the shock wave initiates combustion of the shocked gas. 

E9 The combustion initiated by passage of the shock produces rapid 
thermal expansion which then drives the shock in a self-r e,nforc,n 9 
manner (i.e. the shock wave becomes a detonation wave). 

E10 The detonation wave propagates to the boundary of the cloud of 
explosive gas, after which it degenerates to a (nonreacting) blast 

wave. * 

Ell The blast wave propagates until it strikes the shuttle orbitei . 


5.2 NECESSARY CONTRIBUTING CONDITIONS 


Several of the above events can take place only if certain conditions are 


satisfied. Let t denote the time variable and let t E1 , t E2 , t E3 , ... denote 
the times corresponding to events El, E2, E3. ... . respectively. I propose 

that the necessary contributing conditions corresponding to events fcl t 
listed above are: 


NCC1 (t £ t E2 ) A physical mechanism is avaliable to extinguish the usual 
flame in the exhaust plume. 


NCC2 ( t > t E3 ) The amount of explosive gas and its placement relative 
to the orbiter are such that complete detonation of the material in 
the cloud could produce a damaging overpressure on the orbiter. 

NCC3 ( t 2 t E3 ) Dilution of unburned hydrogen with water (from all 
sources) does not make the mole fraction of h 2 to h 2 o so high or 
low as to exceed the flammability limits for all ternary mixtures of 
Hj, h 2 o , and air. 


*lt is, of course, conceivable that the shuttle orbiter may abut the cloud 
of explosive gas, in which case the detonation wave may strike the orbiter 

directly . 
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NCC4 ( t 2 t E3 ) A mechanism is available to ignite the unburned gas 
in the cloud. 

NCC5 ( t > t zu ) A mechanism is available to induce turbulent mixing 
in the fluid containing the flame front. 

NCC6 ( t > t E7 ) The shock wave generated in E7 is of sufficient strength 
to initiate combustion of the shocked gas. 

NCC7 ( t > t E7 ) The state of the shocked gas remains conducive to com- 
bustion during a long enough time interval for fairly complete com- 
bustion to take place. 

NCC8 ( t > t E10 ) The speed (if any) with which which the center of 
the explosive cloud moves relative to the orbiter is less than the 
speed of the blast wave produced in E10 relative to the cloud. 

NCC9 ( t > t E10 ) The influence of the boundaries is such that the shock 
wave or detonation wave strength does not diminish to safe levels 
before the orbiter is struck. 

Conditions NCC1-NCC9 are 'necessary' in the sense that if any one of them 
is blocked, the sequence of events El -Ell would be interrupted and a catas- 
trophe avoided. The study by Howard ( 1987) 2 cited in the Introduction 
may be interpreted as an effort to evaluate the credibility of NCC2. His 
results indicate that NCC2 is quite credible, i.e. it is likely to be satisfied 
in a typical flight readiness firing or launch abort. The study by Bransford 
S Voth (1987) 4 , which was also cited in the Introduction, may be interpreted 
as an effort to evaluate the credibility of NCC2 and NCC3. Their results 
indicate that neither of these results is easy to discredit, owing primarily 
to the great practical difficulty of estimating the time-dependent spatial dis- 
tributions of the relative concentrations of h 2 , h 2 o. and air in the flame 
trench under realistic conditions. 

The description given in section 3. 2. 1.1 above of how a spherical detonation 
may be initiated suggests that there may be some necessary contributing 
conditions in the above list that are easier to discredit than are the condi- 
tions examined by Howard ( 1987) 2 and Bransford 8 Voth (1987) 4 . Thus, 
for example, the residual jet momentum of the exhaust plume in the flame 
trench may cause enough distortion of the fluid to affect conditions NCC6, 
NCC7, or NCC8. It may well be that the (benign) statistically stationary 
flow in the combustor of a jet engine may be a more realistic model of flow 
in the flame trench than is the model of a transient explosion in a quasi- 
steady body of fluid. 
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6^ JUXTAPOSITION OF REASSURING AND DISQUIETING FACTS 
6.1 REASSURING FACTS 

In assessing the hazard of hydrogen explosions at space shuttle launch pads 
one may identify some facts as reassuring The followinng list is represen- 
tative. 

• If one accepts the idea that the explosivity of a gaseous mixture 
can be quantified (perhaps by the detonation cell size), then by 
any such characterization the sound supression water spray acts 

to make the gases in the exhaust plume of the SSME less explosive. 

• The most sobering experimental observations on deflagration-to-deto- 
nation transition discussed in section 4.3 above were under conditions 
that differ substantially from those that prevail in the flame trench 

at space shuttle launch pads under real operating conditions. In 
particular, flow due to passage of a deflagration through an initially 
stationary well-mixed cloud of explosive gas (as in the most sobering 
experiments) is not the same as flow in the flame trench (which may 
more nearly resemble the flow in a jet engine combustor). 

• The lower density of hydrogen arid water vapor compared to air at 
the same temperature and pressur e provides a mechanism for disper- 
sal of explosive gases, namely bouyant convection. 

• Residual jet momentum in the deflected exhaust plume acts to trans- 
port the explosive gas away from the orbiter . 

• Accidental detonations of gases from h 7 - o 2 rocket engines are rather 
rare events. Thus, Littlefield ( 1987) 3 was able to document only 
six confirmed detonations out of over 16,000 test firings. 

• No explosions have been observed in two aborts and six flight readi- 
ness firings at KSC space shuttle launch pads. 

6.2 DISQUIETING FACTS 

Alongside the above list of reassuring facts, one may list some that are dis- 
quieting. Thus, 

• The amount of unburned hydrogen discharged by the SSMEs after 
the hydrogen burn off ignitors are spent (say 400 pounds) is very 
much greater than the amount needed to produce a blast wave cap- 
able of damaging the orbiter (say 6 pounds if detonated at a distance 
200 feet from the base of the orbiter) [Howard (1987)- ]. 
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• Several experimental studies have led to the conclusion that deflagra- 
tion- to-detonation transition is scale-dependent and that larger scales 
are more conducive to it [Berman (1985) i: ]. 

• Several experimental studies have led to the conclusion that deflagra- 
tion to detonation transition can be triggered by passage of a defla- 
gration across an obstacle and that the threshold size of such an 
obstacle (as a multiple of the detonation cell width A) can be quite 
small [ cf . Sherman, et aJ. (1986) and section 4.1.2 above]. 

• Older experimental studies like the one conducted by the Arthur 

D. Little Company [Anonymous (I960) 25 ] which led to the conclusion 
that unconfined detonations in spherical balloons containing mixtures 
of hydrogen and air could not be produced by weak ignition sources 
(like sparks or squibs) were predicated on the assumption that the 
rate of energy release by the resulting spherical flame is not accel- 
erated by turbulent mixing. Such an experimental test might be 
quite unrepresentative of the situation in a real flame trench, where 
the flow is likely to be turbulent. This admonition holds a- fortiori 
if the trench is loaded with instrument holders and other turbulence- 
producing obstacles. 

• Finally, although accidental detonations of gases discharged by H 7 -o 2 
rocket engines are indeed rare events, they do nevertheless occur. 

7^ RECOMMENDATIONS 

The work that led to this report was guided by the two questions listed in 
section 1.2 above. I propose that the answer to the first question 'Is there 
a set of events common to all credible scenarios leading to a catastrophic 
explosion near the orbiter?' is 'yes'. I have thus proposed a list of such 
a set of events (cf. section 5) and have arranged these events into chrono- 
logical order. To address the second question, '[Are any] of the above 
events precluded by present hardware and operating procedures at KSC?' 

I have proposed a list of 'necessary contributing conditions' all of which must 
be met if the catastrophe scenario cited above is to be credible. The idea 
is that if any one of these necessary conditions is blocked, then the catas- 
trophe is prevented. After reflecting on the credibility of the nine necessary 
conditions listed in section 5, I have not identified a single one which is 
easy to discredit, though some might be less difficult to discredit than 
others. 

In an ideal world, one might aspire to achieve a quantitative scientific under- 
standing of the complete three-dimensional unsteady flow in the flame trench 
under all anticipated operating conditions. The combination of complicated 
boundary geometry, and the multiplicity of physical processes ( e.g. phase 
changes, chemical reactions, heat transfer, etc.) operating at once in the 
flow, however, easily exceeds the capabilities of ali methods for achieving 
such an understanding (computation, analysis, and experiment) with which I 
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am familiar, nor do I believe that such capabilities could be developed in 
the forseeable future, even with the expenditure of copius Resources. \n 
formulating recommendations for further work, therefore one should take 
care to recommend tasks which at least appear to be tractable. For this 
purpose, I propose a third guiding question (which augments the two ques 
tions given in section 1.2 above), namely 

03 What conditions would have to be met before one could, in good con- 

Q science, abandon further work on the unburned hydrogen problem. 

Such a list of conditions might include the following ones. 

Cl One must accept the proposition that there is a set of events common 
to all credible scenarios for a catastrophe and be confident that one 
has identified all of the events in that set. 

C2 One must be confident that one has identified a set of necessary 
contributing conditions to a catastrophe, any one of which, it 
blocked, would forestall the catastrophe. 

C3 One must be confident that there is at least one necessary condition 
that is absolutely precluded by present hardware and operating pro- 
cedures. Here, of course, redundancy is preferable. 

The foregoing discussion leads to the following recommendations. 

• Devote some future effort on the unburned hydrogen problem to fur- 
ther consideration of the list of necessary contributing conditions 
for a catastrophe. 

• After adopting a list of necessary contributing conditions for a catas- 
trophe, devote some future effort to finding and interpreting simple 
ideal models of selected flow details. Such details might include. 

(i) 'flame holding' by stationary obstacles in the flame trench (in 
the manner of 'flame holders' in the combustor of a jet engine J; 

(ii) distortion of an initially spherical flame front by turbulent mining 
and its possible effect on the initiation of spherical detonations; 

(iii) the action of the sound supression water spray on the SSME 
exhaust plume, particularly the manner in which drops of liquid water 
disintegrate and ultimately affect the concentrations of gaseous tf 2 o 
in the region below the spray; (iv) the action of bouyancy and resi- 
dual momentum in the deflected exhaust plume on dispersal of the 
unburned hydrogen. 

Interpretation of simple models of flow details such as those listed above 
might permit one to dismiss one or more of the necessary contributing condi- 
tions for a catastrophe as incredible. Alternatively if such 1 fut V"j ' S V^ '® S 
indicate that certain catastrophe scenarios are indeed credible, then changes 
in hardware and operating procedures inspired by such studies could produce 
a real improvement in the safety of shuttle operations. 
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8. CONCLUSIONS 


• A local fireball whose flame front propagates only subsonicaliy ( i.e. 
a deflagration ) may evolve into one whose flame front propagates 
supersonically (i.e. a detonation) if the rate of release of thermal 
energy is accelerated by turbulent mixing or by any other mechanism 
that rapidly increases the flame area. Such deflagration-to-detona- 
tion transition is possible even in completely unconfined flow. Thus, 
the absence of blasting caps or other high energy sources in the 
KSC flame trench does not, by itself, justify the belief that a hydro- 
gen-air mixture in it could not be detonated by other ignition 
sources. This conclusion follows from the experimental studies cited 
in section 4.3, the theoretical considerations reviewed in chapter 3 
[see especially section 3. 2. 1.1] and it is consistent with the conclu- 
sions of a study by the Arthur D. Little Company on spark ignition 
of hydrogen-air mixtures in spherical balloons [see the penultimate 
'disquieting fact' cited in section 6.2]. 

• Any flammable mixture of gases is detonable in the sense that it is 
possible to produce a one-dimensional detonation wave that propagates 
through it. Indeed, such a detonation wave may always be produced 
by placing the gas in a long tube and driving it with a piston. 
Specifically, if the gas is initially at rest and if the piston undergoes 
a step change in speed from zero to some constant speed u. then 

a detonation wave will always result if u is above a certain threshold 
(whose value depends upon the chemistry of the gas and its initial 
thermodynamic state). In such an experiment, the boundary does 
work on the fluid. ' Detonability limits' can therefore be distinct 
from and narrower than 'flammability limits' only of one adopts a def- 
inition of detonability limits that rules out such working on the fluid 
by the boundaries (at least beyond a certain limited time interval). 
This conclusion follows from the discussion in section 3.1.3 above. 

• Accurate modeling of all the phenomena in the flame trench is not 
possible at present nor is it likely to become possible any time in 
the forseeable future. Even without such accurate modeling, 
however, there are realistic prospects for ruling out accidental hydro- 
gen detonations at space shuttle launch pads. By identifying a list 
of necessary contributing conditions for a catastrophe and by formu- 
lating simple ideal models of selected flow details, there is reason 

to believe that, in time, one may either discredit one or more of the 
necessary contributing conditions (thus certifying the safety of the 
present shuttle system) or identify a change in hardware or opera- 
ting procedures (which would permit such a certification of the modi- 
fied system). Such eforts would seem to be worthwhile. 
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